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Abstract 

A search for CP violation in — ?• (pir^ decays is performed using data collected in 
2011 by the LHCb experiment corresponding to an integrated luminosity of 1.0 fb~^ 
at a centre of mass energy of 7TeV. The CP-violating asymmetry is measured to be 
(-0.04 ± 0.14 ± 0.13)% for candidates with K' K+ mass within 20MeV/c2 of the (/) 
meson mass. A search for a CP- violating asymmetry that varies across the <j) mass 
region of the — )• K^K^tt^ Dalitz plot is also performed, and no evidence for 
CP violation is found. In addition, the CP asymmetry in the Df — t- KgTT~^ decay is 
measured to be (0.61 ± 0.83 ± 0.13)%. 
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1 Introduction 



Cabibbo-suppressed charm decays are the focus of searches for direct CP violation (CPV) 
in the charm sector. In these decays, direct CPV wiU occur if tree and loop (penguin) 
processes interfere with different strong and weak phases. Furthermore, contributions 
from physics beyond the Standard Model may appear in the virtual loops [l]. Evidence 
for direct CPV in charm decays was reported by LHCb and subsequently by CDF using 
the D° — )■ K~ K'^ and — ^ tt^tt^ channels [2|[3]. Although the latest results do not 
confirm the evidence for CPV in the charm sector |4||5], further studies using different 
decay modes remain well motivated. The large branching ratios of — )■ K~ compared 
to — )■ TT"??"*" decays, and of the — )• K~K^'k^ compared to the h'tt^tt^ 
mode, suggest that the contribution of the penguin amplitude may be significant in 
both — 7- K''K~^ and — )■ K^K^n^ decays [g]. The inclusion of charge conjugate 
decays is implied where appropriate throughout this paper. In decays, a non-zero CP 
asymmetry would indicate unambiguously the presence of direct CPV. The — )■ ^vr"*" 
decay is a particularly promising channel for CPV searches due to its large branching ratio 
of (2.65 ± 0.09) X 10~^ I?]. A recent investigation of this decay at the Belle experiment 
yielded a CP-violating charge asymmetry of (+0.51 ± 0.28 ± 0.05)% [s], while BaBar 
measured (-0.3 ± 0.3 ± 0.5)% |9|. 

Searches for CPV in charm decays with the LHCb experiment rely on a good under- 
standing of the charge asymmetries both in D meson production in pp collisions and in the 
detection of the final states. These effects are studied using control decay modes in which 
no CPV is expected, and cancelled by measuring the differences in asymmetries between 
different final states or by comparing measurements made in one area of the Dalitz plot 
relative to another. 

To investigate CPV in the — )■ 07r+ decay, the — )■ K^^tt^ decay with Kg — )■ tt^tt"'" 
is used as a control channel. This decay is itself sensitive to CPV via the interference 
of Cabibbo-favoured and doubly Cabibbo-suppressed amplitudes. However, the CP 



asymmetry in this channel is predicted to be at most 0.01% in the Standard Model 10 
and there is less scope for contributions from non-Standard Model dynamics than in 
the — )■ 071"*" decay as no penguin amplitudes contribute [I]. Therefore CPV in the 
— i- K^TT^ decay is assumed to be negligible. The CP asymmetry in the — > ^vr"*" 
region of the — ^ K^K^n^ Dalitz plot is given by, to first order, 

Acp{D+ ^ 07r+) = A,^{D^ ^ </)7r+) - A,^{D+ ^ i^°7r+) + Acp(ir7^°), (1) 

where the raw charge asymmetry A^aw is defined as 

for yields N£)± of positively- or negatively-charged signal or control-mode candidates. The 
kaon asymmetry Acp{K^ / K^) is the correction for CPV in the neutral kaon system and is 
—0.028% with a systematic uncertainty of 0.028% |11|. To first order, the use of the control 
channel cancels the effects of the production asymmetry of (—0.96 ± 0.26 ± 0.18)% [llj 



1 




m2 (K";i+) [GeV V] 



Figure 1: Variation of the overall phase of the decay amplitude in the (f) mass 
region of the Dalitz plot, from a simulation study based on the CLEO-c amplitude 
model in which the phase is defined relative to that of the i^'*(892)° resonance 14 . To 
calculate Acp\s, the region is divided into rectangular zones as shown, corresponding to 
1.00 < m{K'K+) < 1.02GeV/c2 and 1.02 < m{K'K+) < 1.04GeV/c2 along the y-axis, 
and to m^{K-TT+) < 1.48GeV7c^ and m'^{K-n+) > 1.48GeVyc^ along the x-axis. 



and of any asymmetry associated with the detection of the pion [12] . In the proximity 
of the (j) meson mass of 1019.46 ± 0.02MeV/c2 [t] in the D+ K-K+tt+ Dalitz plot, 
the kaons have almost identical momentum distributions. Therefore the kaon interaction 
asymmetry cancels between the and K~ meson daughters of the resonance. Hence the 
search is restricted to decays with K^K^ invariant masses in the range 1.00 < mK-K+ < 
1.04 GeV/c2. 



A concurrent measurement of the CP asymmetry in the — )■ K^tt^ decay, approxi- 
mated as 



Acp{Dt ^ ir°7r+) = A^^iDj ^ ir°7r+) - A^^iDj ^ 07r+) + Acp(ir7^°), (3) 

is performed using the D'^ — )■ 07r+ decay as a control channel. This decay is also Cabibbo- 
suppressed, with similar contributions from loop amplitudes as the — t- </)7r+ decay, but 
the number of signal candidates is substantially lower. This is partly due to the lower 
production cross-section |13| and partly because only Kg mesons with decay times of 
less than 40 ps are used in this analysis. In Eq. rtsl), the effect of the CPV in the neutral 
kaon system has a sign opposite to that in Eq. ([T| relative to the raw asymmetry in the 
—7- K^TT^ decay because the D'^ decays predominantly to a meson while the 

decays to a K'^. 

Within the Standard Model, CPV in singly Cabibbo-suppressed charm decays with 



contributing tree and penguin amplitudes is expected to be 15 



Im 
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Table 1: Expected mean values of Aqp and Acp\s for different types of CP violation 
introduced into the simulated Dalitz plots, together with the significance with which a 
signal could be observed given estimated overall uncertainties in Acp and Acp\s of 0.2%. 

Type of CPV Mean Aqp (%) Mean Acp\s (%) 

3° in phase -0.01 (0.1a) -1.02 (S.lcr) 

0.8% in ^ amplitude -0.50 (2.5a) -0.02 (0.1a) 

4° in K*{U30f phase 0.52 (2.6a) -0.89 (4.5a) 

4° in ir* (800) phase 0.70 (3.5a) 0.10 (0.5a) 



where i? is a number of order one that depends on hadronic matrix elements, 6s is the 
strong phase difference between tree and penguin amplitudes, and Vij are elements of 
the CKM matrix. In the region of the resonance in the — >■ K~K~^7r~^ Dalitz plot, 
several other amplitudes contribute to the overall matrix element and interfere with 
the meson [9 , 14 . A recent amplitude analysis of this decay mode from the CLEO-c 



collaboration [14] yields a matrix element with a relative strong phase that varies rapidly 
across the region, as shown in Fig. [T| The isobar amplitude model favoured by CLEO-c 
(fit 'B' in Ref. [li]) contains major contributions from the 0, K*{892)°, K^{U30)° and 
-^^^(SOO) resonances. The phase is measured relative to that of the i^'*(892)'' meson. The 
variation in phase means that it is possible that a constant CP-violating asymmetry could 
be cancelled out when the different regions of the resonance are combined to calculate 
Aqp- Hence we define a complementary observable called Acp\s- The area around the 
resonance in the Dalitz plot is split into four rectangular regions A — D defined clockwise 
from the top-left as shown in Fig. [ij The division is chosen to minimise the change in 
phase within each region. A difference between the two diagonals, each made of two 
regions with similar phases, is calculated as 

Acp\s = 2 (^mw + ^mw - ^mw " ^^^w) • i^) 

This observable is not affected by the production asymmetry and is robust against 
systematic biases from the detector. 

To test the hypothesis that Acp\s can sometimes be more sensitive to CP violation 
than Acp, a study is performed using simulated pseudo-experiments in which plausible 



types of CPV are introduced into the CLEO-c amplitude model 14 . The matrix elements 



for and D~ decays are separately modified in a number of ways, as specified in 
Table [l| and events are generated from the resulting probability density functions. In 
each simulated sample, approximately the same number of events as in the dataset are 
produced, and the values of Aqp and Acp\s are compared. The effects of background 
and of the reconstruction and signal selection efficiency variation across the region are 
negligible. 

The level of CPV in the pseudo-experiments is chosen to give an expected result with 
significance of around three Gaussian standard deviations in at least one observable. For 



3 



each type of CPV, twenty Dalitz plots are simulated. The mean values of Aqp and Acp\s 
in these pseudo-experiments are given in Table [Tj together with the significance with which 
these signals could be observed in the dataset under study. The table indicates that some 
types of CPV can be observed more effectively with Acp and others with Acp\ 



s- 



It was found in Ref. 16 that the sensitivity to CPV can vary substantially with 
the details of the amplitude model. Therefore these simple simulations should not be 
treated as accurate predictions, but instead as a guide to the relative sensitivity of the 
two observables. 



2 Detector 



The LHCb detector 17 is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < rj < 5, designed for the study of particles containing 6 or c quarks. The 
detector includes a high precision tracking system consisting of a silicon-strip vertex 
detector (VELO) surrounding the pp interaction region, a large-area silicon-strip detector 
located upstream of a dipole magnet with a bending power of about 4Tm, and three 
stations of silicon-strip detectors and straw drift tubes placed downstream. The combined 
tracking system has momentum resolution Ap/p that varies from 0.4% at 5GeV/c to 0.6% 
at lOOGeV/c, and impact parameter resolution of 20 \mi for tracks with high transverse 
momentum px- Charged hadrons are identified using two ring- imaging Cherenkov detectors. 
Photon, electron and hadron candidates are identified by a calorimeter system consisting of 
scintillating-pad and preshower detectors, an electromagnetic calorimeter and a hadronic 
calorimeter. Muons are identified by a system composed of alternating layers of iron and 
multiwire proportional chambers. The trigger [l8j consists of a hardware stage, based on 
information from the calorimeter and muon systems, an inclusive software stage, which uses 
the tracking system, and a second software stage that exploits the full event information. 



3 Dataset and selection 



The data sample used in this analysis corresponds to an integrated luminosity of 1.0 fb^^ 
of pp collisions at a centre of mass energy of 7TeV, and was collected by the LHCb 
experiment in 2011. The polarity of the LHCb magnet was changed several times during 
the run, and approximately half of the data were taken with each polarity, referred to 
as 'magnet-up' and 'magnet-down' data hereafter. To optimise the event selection and 
estimate background contributions, 12.5 million pp collision events containing — )■ K^ir^, 



— )• TT 7r+ decays and 5 million events containing D 
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with Pythia 6.4 19 with a specific LHCb configuration 
by EvtGen 



— )■ K K^Ti^ decays are simulated 
Hadron decays are described 
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The interaction of the generated particles with the detector and its 



response are implemented using the Geant4 toolkit 22 as described in Ref. 23 



To ensure the dataset is unbiased, the trigger must accept candidates in well-defined 
ways that can be shown to be charge-symmetric. A trigger decision may be based on 
part or all of the D^-j signal candidate, on other particles in the event, or both. For 
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example, signal decays triggered at the hardware level exclusively by the pion from the 



decay are not used, as they are shown in Sect 
charge asymmetries. For an event to be accepted by 



to have large detector-dependent 
jhe hardware trigger, two criteria. 



not mutually exclusive, are used: the decision must be based on one of the daughter 
tracks of the Kg or meson, or on a particle other than the decay products of the -Dj^ 
candidate. In the first case the same track must also activate the inclusive software trigger. 
This software trigger requires that one of the tracks from the signal candidate has 
Pt > 1.7GeV/c and distance of closest approach to the primary vertex (PV) of at least 
0.1mm. The second stage of the software trigger is required to find combinations of three 
tracks that meet the criteria to be signal decays. 

Candidate Dj-j — )■ ^tt"*" decays are reconstructed by combining the tracks from two 
oppositely charged particles that are identified by the RICH detectors as kaons with one 
track identified as a pion. The combined invariant mass of the two kaons is required to 
lie in the range 1.00 < mK-K+ < 1.04GeV/c^. The scalar sum of the of the daughter 
particles must exceed 2.8GeV/c. 

To reconstruct — )■ i^gTr"*" candidates, pairs of oppositely charged particles with a 
pion mass hypothesis are combined to form candidates. Only those with > 700 MeV/ c 
and invariant mass within 35 MeV/ of the world average mass [t] are retained. 
Accepted candidates are then combined with a third charged particle, the bachelor pion, 
to form a D'^^^ candidate. The mass of the meson is constrained to its known value in 
the kinematic fit. All three pion tracks must be detected in the VELO, so only mesons 
with short decay times are used. 

Further requirements are applied in order to reduce background from random track 
combinations and partially reconstructed charm and B decays. Both Kg and -D^^ candi- 
dates are required to have a vertex with acceptable fit quality. Daughters of the and 
Kg mesons must have momentum p > 2 GeV/ c and p^ > 250 MeV/ c. Impact parameter 
requirements are used to ensure that all the daughters of the D'^^^^ candidate do not 
originate at any PV in the event. To remove non-resonant — )■ 7c~TT~^7r~^ candidates, 
the Kg meson decay vertex must be displaced by at least 10 mm in the forward direction 
from the decay vertex of its parent meson. The bachelor pion in both final states 
must have p > 5 GeV/c and p^ > 500MeV/c, must not come from any PV, and must be 
positively identified as a pion rather than as a kaon, electron or muon. In addition, fiducial 
requirements are applied j2j to exclude regions with a large tracking efficiency asymmetry. 
The candidate is required to have 1.5 < < 20.0GeV/c and pseudorapidity r] in the 
range 2.2 < r] < 4.4, to point to a PV (to suppress D from B decays), and to have a decay 
time significantly greater than zero. The proportion of events with more than one Dj-j 
candidate is negligible. 

The invariant mass distributions of selected candidates in the two final states are 
presented in Fig. [2} After applying the selection and trigger requirements, 1,203,930 
Df^^ K^TT^ and 4,704,810 Df^^ 07r± candidates remain in the mass ranges shown in 
the figure. The distribution of decays in the region of the — )■ K^K^n^ Dalitz plot 
is shown in Fig. [3j 
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Figure 2: Invariant mass distribution of selected (a) — t- 07r+, (b) — ?► </)7r^, (c) 
— > f^gTr^ and (d) — t- KgTT~ candidates. The data are represented by symbols 
with error bars. The red dashed curves indicate the signal lineshapes, the green solid 
lines represent the combinatorial background shape, and the green dotted lines represent 
background from mis-reconstructed — )• ^tt+tt" decays in (a) and (b), and — )■ 
KgTT+TT^ or — )• K^K^ decays in (c) and (d). The blue solid lines show the sum of all 
fit components. 



4 Determination of the yields and asymmetries 

For the measurement of Acp-, the signal yields are measured in 12 bins of transverse 
momentum px and pseudorapidity 77, using binned likelihood fits to the distributions of 
the invariant masses m, where m is either m^T^+ or m/^o^+. The values of Acp in each bin 
are calculated and a weighted average over the bins is performed to obtain the final result. 
This procedure is adopted because the distributions of the two decays in and t] differ 
slightly, as shown in Fig. |4| and the production asymmetry may also vary over this 



range 11 . The Pt — V binning therefore reduces a potential source of systematic bias. 
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Figure 3: Observed density of decays in the D 
A-D labelled as described in the text. 
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Figure 4: Distributions of transverse momentum and pseudorapidity rj for (a) — )■ 
K'^TT^ and (b) — )■ ^tt"*" candidates with invariant masses m in the range 1845 < m < 
1895MeV/c^. Candidates that do not fall into the 12 rectangular bins are not used in the 
analysis. 



The shapes of the Z^^^ — )■ Kg7i~^ mass peaks are described by single Cruijff functions 



24 



/(m) oc exp 



— (m — /i)^ 



2cr + (m — jiYaL, 



(6) 



with the peak position defined by the free parameter /x, the width by cr, and the tails by 
ai and a^. The parameter is used for m < /i and for m > /i. In the 07r+ final state, 



Crystal Ball functions 25 are added to the Cruijff functions to account for the tails of 



the mass peaks. The signal lineshapes are tested on simulated data and found to describe 
the data well. The background is fitted with a straight line and an additional Gaussian 
component centred at low mass to account for partially reconstructed — )■ Kg (0)7r+7r° 
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Table 2: Numbers of signal candidates in the four decay modes from the mass fits, with 
statistical uncertainties only. 



Decay mode Signal yield (xlO^) 



1576.9 ± 1.5 
3010.2 ±2.2 
1057.8 ±1.2 
25.6 ±0.2 



decays. This background mostly lies outside the interval in invariant mass that is fitted. 
In the Kg7T~^ case there is also a cross-feed component from the — )■ K^K^ decay mode, 
where the meson is misidentified as a pion. In the fit to data, the cross-feed yield 
and charge asymmetry are allowed to vary but the shape is fixed from the simulation. It 
is modelled by a Crystal Ball function. The yield of cross-feed is found to be small, at 
6014 ± 817 decays, or 0.57% of the D+ yield. 

The fits are performed simultaneously over four subsamples (-Dj^ magnet-up, 
magnet-down, D'^^^ magnet-up, and D^^^ magnet-down data) with the peak positions, 
widths and yields of the and background allowed to vary independently in the four 
subsamples. All other parameters are shared. The peak positions are found to differ 
between charges and magnet polarities by around 0.2MeV/c^. The raw asymmetries are 
then determined from the yields. The fitted yields are given in in Table [2j 

The results are cross-checked with a sideband subtraction procedure under the assump- 
tion of a linear background. The background is sufficiently small relative to the signal in 
the — )■ channel that the charge asymmetry can be calculated by counting and 
D~ candidates in a mass interval defined around the mass of 1845 < m < 1895 MeV/c^. 
Therefore, the yields for Acp\s are evaluated using this simple technique. The resolution 
in the Dalitz plot is improved by constraining the candidate mass to the world average 
value jTj, instead of leaving it as a free parameter. This has a small effect which is assigned 
as a systematic uncertainty in Sect. |5j In the measurement of Acp, the background in the 
— )■ i^gTr"*" channel is larger and therefore the results are taken from fits. 



5 Systematic uncertainties and cross-checks 

The analysis methods are constructed to ensure that systematic biases on the raw charge 
asymmetries cancel in the end result. The dominant systematic uncertainties in both Aqp 
and Acp\s are determined by considering control decay channels in which no asymmetry 
is expected. 

The main systematic uncertainty in Acp results from kinematic differences between 
the 07r+ and Kgir^ final states, which lead to imperfections in the cancellation of detector 
asymmetries between them. Some detector asymmetries arise from small differences in 
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the tracking efficiency or acceptance across tlie bending plane of tlie magnet, i.e. between 
the left and right halves of the detector. The response of the hardware trigger is also 
known to be asymmetric, because it does not take into account which way particles bend 
in the magnetic field when it measures their transverse energy Et or px- In data taken 
with one magnet polarity, a pion from a decay will bend in the opposite direction to 
a pion from a D~ decay, and if one of these pions is bent into an inefficient part of the 
tracking system and is lost while the other is detected, a charge asymmetry will result. The 
same situation could occur if one pion is bent inwards and so does not meet the hadron 
trigger Et threshold while a pion of the opposite charge is bent outwards and therefore 
has enough measured Et to activate the trigger. This cancels to a good approximation 
between the left and right halves of the detector, but any left/right asymmetry in the 
calorimeters or muon stations could result in imperfect cancellation, biasing the charge 
asymmetry. The effect of these asymmetries on this analysis is not eliminated by the 
subtraction of the asymmetries in the two final states in Eq. [l] as the two decays do not 
have identical kinematic properties. Thus, in the data taken with one magnet polarity, 
the charge asymmetry can be affected. However, when the magnet polarity is reversed, 
the bias on the asymmetry changes sign because the particles are deflected in the opposite 
directions. The values of Acp in — )■ ^vr"*" decays are found to differ by (0.81 ± 0.28)% 
between the data taken with magnet polarity up and data taken with polarity down. 
The effect is removed, to a very good approximation, by combining results obtained with 
opposite magnet polarities, Ajjp and A^p, in an average with equal weights, 

^CP - ^ • [I) 

However, non- cancelling effects can bias the measurement and are considered as sources 
of systematic uncertainty. The data triggered by the or meson at the hardware 
level are charge-symmetric to a good approximation, and are assumed to be unbiased. 
However, in data triggered by another particle in the event, the particle that activates 
the trigger may be correlated to the signal decay. For example, a signal decay is often 
accompanied by a meson of the opposite charge. If this meson decays to a charged 
hadron, electron or muon, the daughter particle, which is more likely to have the opposite 
charge to the signal D^, could fire the trigger. The different kinematics and acceptance 
of the signal and control channel studied in this analysis mean that the cancellation 
of charge-asymmetric trigger efficiencies between them may not be complete. To study 
the size of this effect, a sample of approximately 57 million — )■ K^ti^h^ decays is 
selected using the same criteria as those for the signal. The charge asymmetries in the 
differently triggered datasets are given in Table [3} Small but significant discrepancies 
between data from different triggers are observed, indicating that the hardware triggers 
may introduce small biases into the dataset. The large difference between magnet up 
and magnet down data in the sample that is triggered by the muon detectors is due 
to a charge-asymmetric threshold in the detector, but this cancels when the magnet 
polarities are averaged. A systematic uncertainty equal to the maximum deviation from 
the average charge asymmetry of (—2.034 ± 0.014)% in any of the triggers is assigned. 
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Table 3: Raw charge asymmetries, in %, in samples of the — )► K ti^tt^ control decay 
in which a particle not from the signal decay activated various hardware triggers. 

Trigger type Magnet up Magnet down Average Difference 



Hadron 
Muon 
Electron 
Photon 

Overall average 



-2.037 ±0.032 
-2.361 ±0.041 
-2.094 ±0.048 
-1.937 ±0.070 
-2.128 ±0.021 



-1.970 ±0.027 
-1.607 ±0.035 
-2.201 ±0.041 
-2.230 ±0.060 
-1.940 ±0.018 



-2.003 ±0.021 
-1.984 ±0.027 
-2.148 ±0.031 
-2.083 ± 0.046 
-2.034 ±0.014 



-0.068 ± 0.042 
-0.754 ±0.053 
0.106 ±0.063 
0.293 ±0.092 
-0.188 ±0.028 



Table 4: Changes to the final result observed with various alternative kinematic binning 
schemes. The default scheme uses four bins of and three bins of rj. The 

variable is the azimuthal angle around the proton beams. The statistical uncertainties 
are determined by subtracting the uncertainties on the alternative result and the default 
result in quadrature. 



Binning Change in Acp (xlO' 

No binning 8.3 ± 3.7 

12 bins (3 x pt, 4 x D+^^ t]) 0.6 ± 1.7 

48 bins (8 x £)+ p^, 6 x t]) -2.9 ± 1.1 

192 bins (2 x 7r+ p, 8 x 7r+ 0, 4 x p^, 3 x £)+ r]) -2.4 ± 1.1 

180 bins (3 x tt+ p^, 5 x 7r+ r/, 4 x Dj^ p^, 3 x r]) 3.5 ± 2.6 

1440 bins (3 x 7r+ px, 5 x 7r+ r^, 8 x 7r+ 0, 4 x Dj^ p^, 3 x rj) 2.5 ± 1.6 



This occurs in the electron trigger and the difference is 0.114%. The precision with which 
effects cancel between (pn^ and Kg7r~^ final states in the analysis cannot be quantified 
accurately. Therefore the most conservative approach is adopted and no cancellation is 
assumed. 

Residual detector asymmetry differences between the — )■ 07r+ and — > K^7t~^ 
decays due to their different kinematics are studied by applying several different kinematic 
binning schemes to the data. The measured asymmetry is found to be stable with 
variations in the binning, suggesting that the detector asymmetries are small. The results 
are summarised in Table |4} The largest discrepancy in raw asymmetry, as expected, 
results from using no kinematic binning, as this does not account for any variation of the 
production asymmetry across the kinematic region. The next largest difference with 
respect to the baseline binning scheme, of 0.035%, is assigned as a systematic uncertainty 
on the asymmetry due to residual kinematic differences between decay modes. 

The Ac'p\s observable is highly robust against systematic uncertainties. Any effect 
that does not vary across the Dalitz plot will cancel in the subtraction in Eq. [5| and 
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effects that do vary with ir~7r+ or K~ invariant mass across the (f) region will also 
cancel when the regions are combined in the diagonal difference. For example, the 
asymmetry in the interaction of the charged kaons with the detector material would affect 
the asymmetry difference between decays with high and low values of K^tx^ invariant 
mass, which is correlated with the momenta of the kaons. However such effects cancel to 
a good approximation in both observables, as shown below. Only quantities that vary 
between the diagonals of the Dalitz plot region would lead to significant systematic biases 
on Acp\s- To test for the presence of such effects, Acp\s is calculated in the — > ^tt"*" 
control decay, which has similar kinematics to the signal despite the different Dalitz plot 
distributions of the events. The result is (—0.120 ± 0.119)%, which is compatible with 
zero as expected. The statistical uncertainty on this result, added in quadrature to the 
central value, gives a measure of the precision with which detector effects are known to 
cancel. Thus a value of 0.169% is assigned as the main systematic uncertainty in Acp\s- 

The systematic uncertainty due to kaon interaction asymmetries is studied by deter- 
mining the effect on the result of enlarging the size of the K~K^ mass window under 
study. This increases the differences between the momentum spectra of the kaons, which 
increases the effect of the interaction asymmetry because it depends strongly on momentum. 
The consistency of this procedure is checked with simulation studies. The systematic 
uncertainty is found to be 0.031% in Acp for the decay, 0.002% for Acp in the 
decay and 0.009% in Acp\s- 

A systematic uncertainty of 0.056% is associated with the resolution in the Dalitz 
plot variables for Acp\s, due to candidates migrating across the boundaries of the regions 
A — D. This is determined by taking the difference between results before and after the 
mass is constrained to the world average value. This procedure is repeated for Acp-, 
but as expected the systematic uncertainty is much smaller. 

Further small systematic uncertainties arise from the mass fitting, from the calculation 



of the effect of the CPV in the neutral kaon system 11 , from the regeneration of neutral 



kaons via their interaction with detector material 26 , from the choice of fiducial cuts, from 
modelling of the cross-feed in the — )■ K^tt^ decay, and from neglecting the background 
in the calculation of Acp\s- In the simulation, the contribution of D from B decays is 
found to differ between the final states by around 1%, and this leads to another small 
uncertainty since the production asymmetries for B and D decays may differ. 

Other potential sources of systematic uncertainty, such as the difference in selection 
criteria between the two final states, are negligible. The kinematic distributions of daughter 
particles are checked for biases. The variation of the asymmetries with time during the 
data taking period is also considered. The systematic uncertainties are summarised in 
Table [H 

As a further cross check, the difference in raw asymmetry between the -D^ — i- 07r+ and 
— > 7r^7r+7r+ decays is calculated. Since these are both Cabibbo-favoured tree- level 
decays, this quantity is expected to be zero. The — )■ 7r^7r+7r+ decay has reasonably 
similar kinematic properties and a similar yield in our dataset to the — )■ K^n^ decay, 
and the — )■ 07r+ is very similar to the corresponding decay. Thus the kinematic 
differences between the final states in the control decays are similar to those in the 
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Table 5: Systematic uncertainties on the three measurements. The abbreviation n/a is 
used where the systematic effect does not apply. The row labelled 'Backgrounds' represents 
the uncertainty in modelling the cross- feed in Acp and the uncertainty from ignoring the 
background in Acp\s- 



Source 


Acp (D+) \%] 


Acp (D+) \%] 


Arp\^ \%] 


Triggers 


0.114 


0.114 


n/a 


Df control sample size 


n/a 


n/a 


0.169 


Kaon asymmetry 


0.031 


0.002 


0.009 


Binning 


0.035 


0.035 


n/a 


Resolution 


0.007 


0.006 


0.056 


Fitting 


0.033 


0.033 


n/a 


Kaon CP violation 


0.028 


0.028 


n/a 


Fiducial effects 


0.022 


0.022 


n/a 


Backgrounds 


0.008 


n/a 


0.007 


D from B 


0.003 


0.015 


0.003 


Regeneration 


0.010 


0.010 


n/a 


Total 


0.133 


0.130 


0.178 



signal channels. 

The D"^ — ^ 7r^7r+7r+ decay is reconstructed using the same selection as for the signal 
decays. The hardware trigger must be activated by a particle that does not form part of the 
signal decay, or by the vr" meson, or by a random tt^ meson. The resulting sample has a 
large background due to random pions from the primary vertex. To remove this, the regions 
of the -D^ —7- 7r^7r+7r+ Dalitz plot in which one of the pions has a low momentum in the 
rest frame are excluded from the sample by removing the areas of the Dalitz plot below 
the /o(980) resonance. The requirement on t^'tt^ invariant mass m^_^+ > 0.75GeVyc^ is 
applied to both tt'tt^ meson pairs. The mass distribution of the candidates that remain 
is fitted with a Cruijff function in the 12 kinematic bins described in Sect. |4]and the raw 
charge asymmetries in the decay are calculated. 

The weighted average of the raw asymmetry differences in the 12 kinematic bins is 
(0.22 ± 0.12)%. The systematic uncertainty on this is similar to that on the main analysis, 
or 0.13%, so the result differs from zero by 1.3 standard deviations. This discrepancy is 
assumed to be a statistical fiuctuation and no additional uncertainty is assigned. 

Many additional cross-checks and comparisons of the data samples are performed. 
The raw asymmetries are consistent with those observed in the measurements of the 



and production asymmetries 11 12 . The different triggers used in the analysis give 
statistically compatible results. A study of the values of Acp in individual bins gives no 
indication of any dependence on px and rj. The regions A — D used in the calculation of 
Acp\s have fully compatible asymmetries. 
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Figure 5: Invariant mass distribution of selected — t- tt^tt^tt^ decays. The data are 
represented by symbols with error bars. The red dashed peaks indicate the signal decays, 
the green solid lines represent the combinatorial background shape, and the green dotted 
lines represent backgrounds from mis-reconstructed tt^tt^ti^'k^ decays. The blue 

solid line shows the sum of all fit components. 

6 Results and conclusion 

Searches for CP violation in the region of the — )■ K^K^tt^ Dalitz plot and in the 
— 7- K^Ti^ decay mode are performed. The results are 



consistent with existing measurements. The first and third measurements assume negligible 
CP violation effects in the Kgir'^ and — )■ ^Tr"*" control channels, respectively. 

The Acp\s observable is shown to increase the sensitivity of the analysis to certain types 
of CP violation significantly, but there is no evidence for CP violation in either decay. 
This is the most precise analysis of CP violation in the region of the — )■ K^K^-k'^ 
Dalitz plot to date. 

Acknowledgements 

We express our gratitude to our colleagues in the CERN accelerator departments for 
the excellent performance of the LHC. We thank the technical and administrative staff 



Acp{D+ ^ 07r+) 
Acp\s{D+ ^ 0vr+) 
Acp{Dt ^ K'.TT^) 



(-0.04 ±0.14 ±0.13)%, 
(-0.18 ±0.17 ±0.18)%, 
(±0.61 ±0.83 ±0.13)%, 



13 



at the LHCb institutes. We acknowledge support from CERN and from the national 
agencies: CAPES, CNPq, PAPER J and FINER (Brazil); NSFC (China); CNRS/IN2P3 
and Region Auvergne (France); BMBF, DFG, HGF and MPG (Germany); SFI (Ireland); 
INFN (Italy); POM and NWO (The Netherlands); SCSR (Poland); ANCS/IPA (Romania); 
MinES, Rosatom, RPBR and NRC "Kurchatov Institute" (Russia); MinECo, XuntaGal 
and GENCAT (Spain); SNSP and SER (Switzerland); NAS Ukraine (Ukraine); STPC 
(United Kingdom); NSP (USA). We also acknowledge the support received from the ERG 
under FP7. The Tierl computing centres are supported by IN2P3 (Prance), KIT and 
BMBP (Germany), INPN (Italy), NWO and SURP (The Netherlands), PIC (Spain), 
GridPP (United Kingdom). We are thankful for the computing resources put at our 
disposal by Yandex LLC (Russia), as well as to the communities behind the multiple open 
source software packages that we depend on. 



References 

Y. Grossman, A. L. Kagan, and Y. Nir, New physics and CP violation in singly Cabibbo 



[1] 
[2] 
[3 

[4 

[5] 
[6] 
[7 



[9 



suppressed D decays, Phys. Rev. D75 (2007) 036008', arXiv : hep-ph/0609 178 



LHCb collaboration, R. Aaij et ai, Evidence for CP violation in time-integrated 



D° h-h+ decay rates, |Phys. Rev. Lett. 108 (20 12) 111602 



arXiv: 1112.0938 



GDP collaboration, T. Aaltonen et ai, Measurement of the difference of CP -violating 



asymmetries in K^K and — )• tt+tt decays at CDF, Phys. Rev. Lett. 109 



1(2012) 11180H|arXiv: 1207 .2158 



LHCb collaboration, R. Aaij et ai. Search for direct CP violation in h h^ 

modes using semileptonic B decays, arXiv: 1303.2614, submitted to Physics Letters 
B. 

LHCb collaboration, R. Aaij et ai, A search for time-integrated CP violation m 



D° ^ K-K+ and n-n+ decays, |LHCb-CQNP-2013-003 

J. Brod, Y. Grossman, A. L. Kagan, and J. Zupan, A consistent picture for large 
penguins m D° n+n-, K+R-, JHEP 1210 (2012) 161 



arXiv: 1203.6659, 



Particle Data Group, J. Beringer et ai. Review of particle physics^ Phys. Rev. D86 
p012) 010001, 

Belle collaboration, M. Staric et ai, Search for CP Violation in Meson Decays to 



Phys. Rev. Lett. 108 (2012) 071801, ,arXiv: 1110.0694, 



BaBar collaboration, J. Lees et ai. Search for direct CP violation in singly- Cabibbo sup- 
pressed K+K'Tc^ decays, |Phys. Rev. D 87 (2012) 052010t |arXiv: 1212. 1856] 



14 



[10] I. I. Bigi and H. Yamamoto, Interference between Cabibbo allowed and doubly for- 
bidden transitions m K^,K^ + tc decays, Phys. Lett. B349 (1995) 363| 
|arXiv : hep-ph/9502238| 

[11] LHCb collaboration, R. Aaij et ai, Measurement of the production asymmetry in 
7TeV pp collisions, Phys. Lett. B718 (2013) 902 , arXiv: 1210.4 112. 

[12] LHCb collaboration, R. Aaij et ai, Measurement of the Df- D~ production asymmetry 
m 7 TeV pp collisions, Fhys. Lett. B713 (2012) 186, arXiv: 1205 . 0897 . 

[13] LHCb collaboration, R. Aaij et ai. Prompt charm production in pp collisions at 
sqrt(s)=7 TeV, ,Nuclear Physics B871 (2013) 1, |arXiv: 1302 .28641 

[14] CLEG collaboration, P. Rubin et ai. Search for CP violation in the Dalitz plot analysis 



of D+ K-K+7r+, |Phys. Rev. D78 (2008) 072003^ |arXiv:0807.4545[ 



[15] D. Atwood and A. Soni, Searching for the origin of CP violation in Cabibbo suppressed 
D-meson decays, [arXiv: 1211 . 1026| 

[16] LHCb collaboration, R. Aaij et ai. Search for CP violation m — )■ K^t[^ 



decays, [Phys. Rev. D84 (2011) 112008| arXiv: 1110.3970] 



[17] LHCb collaboration, A. A. Alves Jr. et al. The LHCb detector at the LHC, |JINST3 
[(2008) S08005. 



[18] R. Aaij et ai. The LHCb trigger and its performance, arXiv: 1211 .3055', submitted 
to JINST. 



[19] T. Sjdstrand, S. Mrenna, and P. Skands, PYTHIA 6.4 physics and manual, JHEP 05 
,(2006) 026, ,arXiv : hep-ph/0603175, 

[20] L Belyaev et ai. Handling of the generation of primary events in Gauss, the LHCb 
simulation framework. Nuclear Science Symposium Conference Record (NSS/MIC)| 
[IEEE (2010) 1155, 



[21] D. J. Lange, The EvtCen particle decay simulation package, Nucl. lustrum. Meth 



A462 (2001) 152 



[22] GEANT4 collaboration, J. Allison et ai, Geant4 developments and applications, 
JEEE Trans. Nucl. Sci. 53 (2006) 270; GEANT4 collaboration, S. Agostinelh et al, 
GEANT4: A simulation toolkit, |Nucl. lustrum. Meth. A506 (2003) 250[ 

[23] M. Clemencic et ai, The LHCb simulation application, Gauss.' design, evolution and 
experience, [j. of Phys: Conf. Ser. 331 (2011) 032023[ 

[24] BaBar collaboration, P. del Amo Sanchez et ai. Study of B ^ X7 decays and 
determination of \Vtd/Vts\, W^YS- Rev. D82~(2010) 05ir01[ [arXiv: 1005 .40871 



15 



[25] T. Skwarnicki, A study of the radiative cascade transitions between the Upsilon-prime 
and Upsilon resonances, PhD thesis, Institute of Nuclear Physics, Krakow, 1986, 
IDESY-F31-86-n2l 

[26] B. Ko, E. Won, B. Golob, and P. Pakhlov, Effect of nuclear interactions of 
neutral kaons on CP asymmetry measurements, ,Phys. Rev. D84 (2011) 11150l[ 
larXiv: 1006 .19381 



16 



